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It is shown that the band structure of single-wall semiconducting carbon nanotubes (CNT) is 
analogous to relativistic description of electrons in vacuum, with the maximum velocity u= 10 8 cm/s 
replacing the light velocity. One-dimensional semirelativistic kinematics and dynamics of electrons 
in CNT is formulated. Two-band k-p Hamiltonian is employed to demonstrate that electrons in CNT 
experience a Zitterbewegung (trembling motion) in absence of external fields. This Zitterbewegung 
should be observable much more easily in CNT than its analogue for free relativistic electrons in 



PACS numbers: 73.22.-f 73.63.Fg 78.67.Ch 

Since the first reported observations of carbon nan- 
otubes [1,2] these unique one-dimensional nanostructures 
were subject of very intensive research, both because 
of their remarkable properties as well as their potential 
use in nanometer-sized electronics (see Refs. 3 and 4). 
Among other particularities CNT have interesting energy 
band structure and it is this aspect that is of our concern 
here. We will be interested in the simplest single- wall 
semiconducting CNT. Such tubes are obtained from a 
slice of graphene wrapped into a seamless cylinder, so 
that the ID band structure of CNT can be constructed 
by using the 2D band structure of graphene. The peri- 
odic boundary condition around the tube circumference 
causes quantization of the transverse wave vector com- 
ponent k x . The purpose of our work is to predict new 
properties of electrons and holes in CNT in the classical 
and quantum domains. To this end we use a similarity 
of the band structure of CNT to the relativistic descrip- 
tion of free electrons in vacuum. In particular, we predict 
that the " semirelativistic" band structure of CNT should 
result in Zitterbewegung (trembling motion) of electrons 
in absence of external fields. Similar phenomenon was 
predicted for relativistic electrons in vacuum but never 
observed. Thus, an observation of this effect in CNT 
would be of great importance not only for the solid state 
physics but also for the high energy physics. 

In the following we use the k-p band structure at the 
K point of the Brillouin zone [5,6]. The initial 2x2 
Hamiltonian is written in the form 



H = a 



a n -ip 
a n + ip 



(1) 



where a is a coefficient, p is the operator of pseudomo- 
mentum in the y direction and a n is given by the quan- 
tization of the wave vector k x . There is a n — hk x (n) — 
h{2ir/L){n - i//3) for n = 0,±1,±2.... Here L is the 
length of circumference. Semiconducting CNT of our 
interest are obtained for v = ±1. In absence of ex- 
ternal fields the resulting energy is £ = ±E(p), where 
E(p) — a{a 2 t -l-p 2 ) 1 / 2 . The upper sign is for the con- 



duction and the lower for the valence band. The above 
relation is analogous to the dispersion E{p) for free rela- 
tivistic electrons in vacuum. We write the energy in the 
following equivalent form 
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where e g = 2aa n is the energy gap and to*, = a n /a is the 
effective mass at the band edge, related to band's curva- 
ture for small p values. Both e g and have different 
values for different subbands n. Equation (2) has the rel- 
ativistic form with the correspondence: e g —> 2m* ) c 2 and 
to*, -> to [7] ■ 

We first investigate some consequences of the one- 
dimensional dispersion (2). The electron velocity is 
v = clE I 'dp = a 2 p/E. For large momenta the velocity 
reaches a saturation value u — a — (e g /2ml) 1 ^ 2 , the 
same for all subbands. The maximum velocity u plays 
for electrons in CNT the role of the light velocity c for 
relativistic electrons in vacuum. Using u we can write 
the energy in another equivalent form 



E(p) = [(toqU 



2\2 



2 2ll/2 
U p \ 1 , 



(3) 



which is directly reminiscent of the relativistic E(j>) rela- 
tion. Now we define an energy-dependent effective mass 
to* relating velocity to momentum 



m v = p. 

We calculate to* = p/v — E/u 2 . This gives 

TP * 2 

E = to u , 



(4) 



(5) 



which is in one-to-one correspondence to the Einstein re- 



lation E = 



(the maximum velocity u replacing c). 



For p — Eqs. (3) and (5) reduce to E r = m^u 2 , 
which corresponds to the formula for the rest energy. 
One can also express the mass to* by the velocity. Be- 
ginning with the relation p 2 u 2 = v 2 E 2 /u 2 , employing 
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E 2 = (m u 2 ) 2 +u 2 p 2 and solving for momentum, we get 



p = m 7U, 



(6) 



where 7 = (1 — v 2 /u 2 )^ 1 / 2 . Using the above definition 
p = m*v we have m = to*,7, which has the familiar rel- 
ativists form (with u replacing c). Now Eq.(5) reads 
E = rriQ-fu 2 . 

Next we assume, in analogy with the special relativity, 
that dp/dt — F, where F is the force. One can now 
define an another effective mass M* , relating force to 
acceleration 



M*a = F. 



(7) 



Since a = dv/dt = (dv/ dp) (dp/dt) = (d 2 E/dp 2 )F, we 
obtain 1/M* = d 2 E/dp 2 . Using the above dispersion 
E(p) we get 



(8) 



which again has the corresponding relation in special rel- 
ativity if the acceleration is parallel to the force. 

Estimating the introduced quantities, we use the value 
of ah = 6.46 eVA [6] . This gives for the maximum ve- 
locity u = a = 0.98 x 10 8 cm/s. The lowest energy gap 
is e g (0) = 2aao, where ao = h2ir/3L. For the circum- 
ference L — 60 A we get e 9 (0) = 0.45eV. The effective 
mass mj/mo = ao/amo = 0.04f for the same condi- 
tions. The quoted parameters are close to those of the 
the typical narrow gap semiconductor InAs, but CNT of 
higher diameter have smaller e g and mj. We emphasize 
that, while for a nonparabolic energy band one can le- 
gitimately define both energy dependent masses to* and 
M* , it is the "momentum" mass to* that is a more useful 
quantity. This mass is employed in the transport the- 
ory since the current is related to velocity, not to ac- 
celeration. It is to* that enters into the equivalence of 
the mass and the energy in Eq.(5). Finally, to* is re- 
lated to the density of electron states since the latter is 
p(E) = (l/nh)(dp/dE) = (1/Trh)(m* /p). 

Next we consider the quantum effects in CNT related 
to the Hamiltonian (1). To this purpose we introduce an 
important quantity [8] 



\z = 



h 



h_ 

a n 



(9) 



which we call the length of Zitterbewegung (see below). 
It corresponds to the Compton wavelength A c for elec- 
trons in vacuum and it plays for the semirclativistic band 
structure (1) the role that A c does for the Dirac equation. 
However, Xz is several orders of magnitude larger than 
A c . Using to*, = 0.041to and u = 0.98 x 10 8 cm/s we 
calculate Xz = 28. 6A. 



Now let us consider the operator of electron velocity 
v = dH/dp. A simple calculation shows that the eigen- 
values of v are, paradoxically, ±u. This differs drasti- 
cally from the classical velocity calculated above. To 
clear the paradox we calculate v(t). Using Eq. (I) we 
get vH + Hv = 2a 2 p = 2u 2 p. Hence the time derivative 
of v is 



dv i „ , „ 
* = H 2UP 



—vH. 

h 



(10) 

This represents a simple differential equation for v. Its 
solution is 



H 



u 2 p^ , Ht . 
-£)exp(-2i—), (If) 



where 1/H = E~ 2 H. Thus the quantum velocity differs 
from the classical one by the term that oscillates in time. 
Equation (11) can be integrated with respect to time to 
give the position operator y in the Hciscnberg picture 

u 2 p ihu ' r ,—2iHt, 
y(t) = y(0) + -^t+^A [exp(^—) - 1], (12) 

where Aq = (vq/u) — (up/H). The first two terms of 
Eq. (12) represent the classical electron motion. The 
third term describes time dependent oscillations with a 
frequency u>z = £g/fr- Since Ao ~ 1 the amplitude of 
oscillations is 2hu/2H w S/to-qU = Xz- In the relativistic 
quantum mechanics (RQM) of free electrons the anal- 
ogous oscillations were devised by Schroedinger [9,10], 
who called them Zitterbewegung (ZB). This explains the 
name given above to Xz- We note that the phenomenon 
of ZB goes beyond Newton's first law since we have a non- 
constant velocity without a force. In RQM it is demon- 
strated that the ZB is a result of interference between 
states of positive and negative electron energies [11]. The 
ZB described above for CNT has the frequency of about 
10 7 times lower and the amplitude of about 10 4 — 10 5 
times higher than the corresponding values for relativis- 
tic electrons in vacuum. We come back to the problem 
of ZB below. 

In order to investigate the case of a definite sign of en- 
ergy we complete the Hamiltonian (1) by the potential 
V(y) on the diagonal and introduce magnetic field by 
adding the vector potential to the momentum operator. 
We then try to separate the 2x2 eigenenergy problem 
with the Hamiltonian (1) into two independent problems 
for each band. This can be easily done in the semiclas- 
sical approximation neglecting the noncommutativity of 
p with V. One obtains then, as expected, the one-band 
effective mass approximations for the conduction and va- 
lence bands 



€ P 2 

H = ±-^± — 

2 2m* 



+ V, 



(13) 



where P = p y 



eA y is the canonical momentum. 
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In absence of external potentials the two-component 
wave functions resulting from the Hamiltonian (1) can 
be transformed exactly into one-component functions for 
positive (or negative) electron energies. This is achieved 
by applying the following unitary transformation S 

where = ^ . It can be verified that SS + = 1 
and that the transformed Hamiltonian is 

Hj, = SHS+ = E p $. (15) 

Thus the positive and negative energies are separated. 
One can transform any wave function ip(y') (not only the 
cigenfunctions of the Hamiltonian (1)) from the original 
two-component representation to the one-component <j> 
representation by the following integral transformation 
(cf. Rcfs. 12,13,14) 

4>'±(y) = J K ± {y,y')i>(y')dy' , (16) 

where 

K ± (y,y') = ±-(1±0) J {I + ^)e*'(»-»V. (17) 

The subscripts ± correspond to the one-component func- 
tions related to positive or negative energy, respectively. 
The factors (1 ± ($) guarantee this property. The kernels 
K±(y,y') are not point transformations. Suppose we are 
interested in the eigenfunction of the electron position y. 
It is convenient to take in the initial representation the 
unit vector tp(y') = Q ) )S(y' — y ). It then follows from Eq. 
(16) that the transformed one-component functions are 
K±{y-yo). The integrals of K± {y-y ) over Y — y — y 
are unity. To get an idea of the extension of K±(y — y ) 
we calculate their second moment. After some manipu- 
lation we obtain 

J Y 2 K ± (Y)dY = \x\. (18) 

Thus the extension of the transformed eigenfunctions of 
position is \y — ya\ — \z/V%- In the transformed states 
there is no ZB since the transformation (14) [ or (16) 
] eliminates the negative (positive) energy components 
of the wave functions. Following the interpretation es- 
tablished in RQM we are confronted with the choice be- 
tween a point-like electron described by a two-component 
function (four components including spin) which experi- 
ences the ZB with the amplitude of Xz, and an elec- 
tron described by a one-component function (two com- 
ponents including spin) for cither positive (or negative) 
electron energy which does not experience the ZB but 



is "smeared" in the y direction to an object of the size 
Xz- One can say that in the one-component represen- 
tation the trembling motion is averaged into smearing. 
As follows from the above estimation of Xz, the ampli- 
tude of ZB or, alternatively, the smearing of electrons in 
CNT is quite large and it should be observable directly 
or indirectly. 

In a recent paper [15] Zitterbewegung-type oscillations 
were proposed using the Hamiltonian of spin splitting due 
to structure inversion asymmetry (SIA) or bulk inversion 
asymmetry (BIA) of the system [16,17,18]. The Hamil- 
tonian for BIA has the form Hd = (f3/h)(a y p y — a x p x ), 
where a x and a y are the Pauli spin matrices. Inter- 
estingly, the Hamiltonian (1) we use for CNT has the 
identical form, with the minus sign in Hp replaced by 
the plus sign in Eq. (1). This similarity is formal be- 
cause the Hamiltonian (1) is not related to spin. The 
description of Ref. [15] begins with a 2D case but fi- 
nally treats a quantum wire, which is a system almost 
identical with CNT. Considering the motion of a wave 
packet it shows that the ZB occurs in the direction per- 
pendicular to the packet's group velocity. Its frequency 
is u) = AE/h, where AE is the energy splitting, and its 
amplitude is inversely proportional to the packet's wave 
vector. Our Eq. (12) describes the same result for the 
frequency, namely u>z = e g/f i - As to the amplitude, our 
result is formally also the same since the ZB in the y di- 
rection has the amplitude Xz — h/a n — l/k x (n), where 
k x {n) is the wave vector in the x direction. However, 
the latter similarity is only formal because the quantized 
value of k x corresponds to the standing wave, that is to 
both +k x (n) and — k x (n). Thus one should rather have 
the packet propagating in the free direction y with the 
wave vector k y o- The ZB is then parallel to the x direc- 
tion, i.e. it occurs on the tubes's circumference, similarly 
to the wire's finite width considered in [15]. 

As compared to the oscillatory electron motion pro- 
posed in [15], the ZB in CNT described above has two 
advantages. First, in the configuration discussed in Ref. 
[15] the growth direction creating the initial 2D system is 
parallel to z. It is known that, for the Bychkov-Rashba 
Hamiltonian related to SIA, the electron spin is quan- 
tized in the plane transverse to the growth direction, 
i.e. in the x — y plane (see Ref. [18]). However, the 
calculation [15] requires that the spin be directed along 
the z direction which contradicts the above quantization 
condition. In our case there is no requirement concern- 
ing the spin direction since the mechanism of ZB is not 
related to spin. The second advantage is more funda- 
mental. The ZB proposed above for electrons in CNT is 
the "true" Zitterbewegung in a sense that it corresponds 
to the ZB for free relativistic electrons in vacuum, while 
the oscillatory motion proposed in [15] has no such corre- 
spondence. Thus the observation of ZB in CNT would be 
of great value also for the relativistic quantum mechan- 
ics. There exist in the literature contradictory statements 
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concerning the observability of ZB in vacuum (see e.g. 
Refs. 10,19,20 ), which makes its analogue in solids and 
molecules even more interesting. It appears that the ZB 
in CNT can be observed experimentally using the high- 
resolution scanning-probe microscopy [21,22]. 

In summary, using the analogy between the band struc- 
ture of single- wall semiconducting carbon nanotubes and 
the description of relativistic electrons in vacuum we for- 
mulate the one-dimensional semirelativistic kinematics 
and dynamics for charge carriers in CNT. We also con- 
sider the quantum domain demonstrating that electrons 
in CNT experience the Zitterbewegung (trembling mo- 
tion) even in absence of external fields. If the electrons 
are described by an effective one-band Hamiltonian for 
positive (or negative) energy, there is no ZB but the elec- 
trons should be treated as objects extended along the 
direction of CNT axis. It is emphasized that the analo- 
gous effects have been predicted for relativistic electrons 
in vacuum but they are much more difficult to observe 
than in CNT. 

I am pleased to thank Dr. T.M. Rusin and Dr. P. 
Pfeffer for elucidating discussions. This work was sup- 
ported in part by The Polish Ministry of Sciences, Grant 
No PBZ-MIN-008/PO3/2003. 



[1] S. Iijima, Nature (London) 354, 56 (1991). 

[2] Z.Ya. Kosakovskaya, L.A. Chcrnozatonskii, and E.A. Fe- 

dorov, JETP Lett. 56, 26 (1992).[Pis'ma Zh. Eksp. Teor. 

Fiz. 56, 26 (1992)]. 
[3] Carbon Nanotubes, edited by M.S. Dresselhaus, G. Dres- 

selhaus, and Ph. Avouris, (Springer, Berlin, 2001). 
[4] R. Saito, G. Dresselhaus, and M.S. Dresselhaus, Physical 



Properties of Carbon Nanotubes (Imperial College Press, 
London, 1999). 

[5] J.C. Slonczewski and R.R. Weiss, Phys. Rev. 109, 272 
(1958). 

[6] H. Ajiki and T.Ando, J. Phys. Soc. Jpn 62, 1255 (1993); 

H. Ajiki and T. Ando, ibid 64, 4382 (1995); T. Ando 

and T. Nakanishi, ibid 67, 1704 (1998); T. Ando and H. 

Suzuura, ibid 71, 2753 (2002). 
[7] W. Zawadzki, in High Magnetic Fields in the Physics of 

Semiconductors II, edited by G. Landwehr and W. Ossau 

(World Scientific , Singapore, 1997), p. 755. 
[8] W. Zawadzki, Phys. Rev. B 72, 085217 (2005). 
[9] E. Schroedinger, Sitzungsber. Preuss. Akad. Wiss. Phys. 

Math. Kl. 24, 418 (1930); Schroedinger's calculation is 

reproduced in Ref. 10. 
[10] A.O. Barut and A.J. Bracken, Phys. Rev. D 23, 2454 

(1981). 

[11] W. Greiner, Relativistic Quantum Mechanics (Springer, 
Berlin, 1994). 

[12] L.L. Foldy and S.A. Wouthuysen, Phys. Rev. 78, 29 
(1950). The RHS of Eq. (21) in this paper should be 
completed by a factor of 1/2. 

[13] M.E. Rose, Relativistic Electron Theory (Wiley, New 
York, 1961). 

[14] J. A. Eisele, Modern Quantum Mechanics with Applica- 
tions to Elementary Particl Physics (Wiley, New York). 

[15] J. Schliemann, D. Loss, and R.M. Westervelt, Phys. Rev. 
Lett. 94, 206801 (2005). 

[16] Y.A. Bychkov and E.I. Rashba, J. Phys. C. 17, 6039 
(1984). 

[17] G. Dresselhaus, Phys. Rev. 100, 580 (1955). 
[18] W. Zawadzki and P. Pfeffer, Semicond. Sci. Technol. 19, 
Rl (2004). 

[19] K. Huang, Am. J. Phys. 20, 479 (1952). 

[20] J.A. Lock, Am. J. Phys. 47, 797 (1979). 

[21] M.A. Topinka et al. Science 289, 2323 (2000). 

[22] B. J. LeRoy, J. Phys. Condens. Matter 15, R1835 (2003). 



